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ABSTRACT 

Ground  truth  measurements  of  the  atmospheric  aerosol  optical 
thickness  and  radiance  data  for  the  GOES-1  and  SMS- 2  satellites,  show  that 
the  satellite  sensors  were  stable  over  a  period  of  eighteen  months,  and 
that  there  is  probably  a  difference  in  the  radiometric  calibrations  of  the 
sensors.  Theoretical  calculations  show  that  errors,  due  to  variations  in 
the  aerosol  properties,  in  inferring  the  aerosol  content  from  a  satellite 
radiance  measurement  are  similar  for  the  nadir  viewing  sensors  (e.g., 
Landsat  MSS)  and  for  the  scanning  sensors  (e.g.,  NOAA-6  AVHRR). 
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1 .  INTRODUCTION 


The  investigation  of  the  satellite  technique^  to  measure  tropo¬ 
spheric  aerosols  has  continued  both  experimentally  and  theoretically. 
Additional  ground  truth  measurements  and  satellite  measurements  have  been 
obtained  for  the  60ES-1  and  SMS-2  satellites,  and  compared  with  earlier 
results.  A  discussion  of  the  relative  sensitivities  of  the  different 
satellite  sensors  for  monitoring  aerosol  changes  is  given.  The  results  of 
calculations  to  investigate  the  errors  in  the  technique  due  to  variations 
in  the  refractive  index,  the  size  distribution,  and  the  vertical  distribu¬ 
tion  of  the  aerosols  are  discussed. 
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2.  GROUND  TRUTH  MEASUREMENTS 


The  acquisition  of  ground  truth  measurements  of  the  aerosol  optical 
thickness  at  the  time  of  satellite  overpasses  has  continued  when  possible 
during  other  Navy  programs.  Measurements  have  been  made  by  NRL  personnel  at 
San  Nicolas  Island,  California;  at  Panama  City,  Florida;  at  Wallops  Island, 
Virginia;  in  the  North  Sea  (Marsen  experiment);  and  on  the  USNS  Hayes  off  the 
Virginia  coast;  in  addition,  measurements  have  been  made  by  SAI  at  San  Diego. 

In  order  to  obtain  reliable  measurements  of  aerosol  optical  thick¬ 
ness  it  is  important  to  have  a  stable  well-calibrated  instrument.  A  recent 

( p ) 

analysis v  '  of  seven  NOAA  photometers  showed  that  over  a  period  of  about  one 
year  the  instrument  calibration  constant  changed  by  as  much  as  16.5%,  and 
that  much  of  the  change  is  probably  due  to  accumulation  of  dirt  on  the 
filters.  The  analysis  suggested  that  if  the  instruments  are  kept  in  a  clean 
dust- free  environment  then  the  uncertainty  in  the  instrument  calibrations 
could  be  kept  within  acceptable  limits  of  a  few  percent. 

In  light  of  the  findings  for  these  NOAA  sunphotometers  (manufactured 
by  Eppley  Laboratories  under  NOAA  specifications)  it  was  decided  to  investi¬ 
gate  the  stability  of  the  Volz  sunphotometers  used  by  NRL  in  this  program.  A 
check  on  the  stability  of  a  sunphotometer  is  best  accomplished  by  periodically 
making  careful  Langley  plots,  or  by  periodically  comparing  it  with  a  reference 
sunphotometer  which  is  regularly  calibrated.  In  the  absence  of  such  checks, 
the  use  of  two  sunphotometers,  as  done  by  NRL,  can  be  very  useful.  If  a  meas¬ 
urement  is  repeated  with  the  second  sunphotometer  immediately  after  making 
one  with  the  first  sunphotometer,  a  plot  of  the  readings  can  be  made,  as  shown 
in  Fig.  1  for  the  NRL  sunphotometers  No.  448  and  No.  457.  This  plot  shows 
that  in  the  period  June  1979  to  April  1980,  the  response  of  the  500  nm  channels 
have  not  changed  with  respect  to  each  other,  indicating  that  the  calibration 
constants  are  probably  unchanged  in  that  period  of  time  (it  is  possible,  but 
unlikely  that  each  has  changed  by  the  same  percentage).  A  similar  plot  for  the 
440  nm  channels  is  given  in  Fig.  2,  and  for  these  channels,  a  changa  in  one  of 
the  instruments  may  have  occurred  between  June  1979  and  April  1980,  although 
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#448  Meter  Deflection 


Figure  2.  Comparison  of  NRL  Sunphotometer  (440  nm) 


three  of  the  earlier  points  are  so  close  to  the  later  points,  a  change  is 
not  obvious.  Thus,  in  analyzing  the  measurements,  we  will  assume  that  no 
significant  change  in  the  calibration  constants  has  occurred  in  this  time 
period. 

Further  evidence  of  the  stability  of  these  sunphotometers  is  pro¬ 
vided  by  data  taken  at  Wallops  Island  on  25  June  1979,  when  a  few  measure¬ 
ments  were  obtained  over  a  range  of  air  masses  sufficient  to  make  rough 
Langley  plots  for  each  instrument,  as  shown  in  Fig.  3.  These  data  are  too 
few,  (and  the  atmospheric  optical  thickness  does  not  seem  constant  over  the 
whole  data-taking  period)  to  make  a  reliable  calibration.  However,  the 
calibration  constants  (I0)  (see  Table  1)  which  would  be  inferred  from  the 
intercepts  on  the  ordinate  are  very  close  to  those  originally  given  with 
instruments,  except  for  the  500  nm  channel  of  No.  448.  Thus,  except  tor 
this  channel,  we  assume  that  the  original  I  values  are  still  valid  through 
April  1980.  For  the  500  nm  channel  of  No.  448  we  assume  I  =  51.5  (deter¬ 
mined  from  Fig.  1,  assuming  I  =  50.0  for  No.  457).  All  the  data  taken  in 
1975  and  1980  indicate  that  the  IQ  for  No.  448  is  larger  than  that  for  No. 
457,  in  contrast  to  the  original  values,  suggesting  either  an  incorrect 
original  value  for  No.  448,  or  a  change  since  it  was  obtained. 
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Air  Mass  (m) 


Figure  3.  Langley  Plots  for  NRL  Sunphotometers 


Table  1.  Sunphotometer  Calibration  Constants 


500  nm 


Given  (1978) 


Measured  (1979) 


Given  (1978) 


440  nm 

Measured  (1979) 


No.  448  47.8  51.1  57.5  58.4 

No.  457  50.0  49.6  39.5  40.8 
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3.  SATELLITE  MEASUREMENTS 

The  original  GOES-1  results^  obtained  during  the  1977  EOMET 
cruise  across  the  North  Atlantic  have  been  supplemented  by  two  more  data 
points  obtained  just  off  the  coast  from  Panama  City,  Florida  in  December 

1978.  The  results  given  in  Fig.  4  show  very  good  agreement  with  the 
Atlantic  data  obtained  nineteen  months  earlier,  suggesting  that  the  sensi¬ 
tivity  of  the  GOES-1  sensor  is  quite  stable  over  long  periods. 

Data  for  the  SMS-2,  the  same  type  of  satellite  as  GOES-1 ,  but 
positioned  over  the  Pacific  Ocean,  were  obtained  for  five  days  at  San 
Nicholas  Island  in  May  1978,  and  for  five  days  at  San  Diego  in  November 

1979.  At  San  Nicolas  Island,  the  aerosol  content  was  effectively  the  same 
each  day,  and  the  SMS-2  measured  effectively  the  same  radiance  each  time, 
thus  verifying  the  repeatability  of  the  technique.  The  data  for  San  Diego 
are  also  given  in  Fig.  4  and  show  that  the  slope  of  the  radiance  vs  aerosol 
content  relationship  for  SMS-2  is  the  same  as  that  found  for  GOES  during 
the  EOMET  cruise.  The  difference  between  the  two  lines  could  be  due  to  a 
difference  in  the  aerosol  properties  between  the  Atlantic  and  Pacific  Oceans, 
but,  as  discussed  previously,^  it  is  most  likely  due  to  a  difference  in 
the  radiometric  calibrations  of  the  satellite  sensors. 

It  is  noted  in  Fig.  4  that  one  of  the  San  Diego  data  points  is 
significantly  higher  than  expected.  However,  if  the  sensor  digital  count  at 
this  time  had  been  just  one  count  less,  then  the  radiance  would  be  1.07  rather 
than  1.42  mW/cm  /ym/Sr,  bringing  it  almost  to  the  expected  value.  This  large 
change  in  radiance  for  a  change  of  one  digital  count,  which  is  certainly  a 
possible  error  in  the  GOES/SMS  system,  clearly  illustrates  the  insensitivity 
of  this  sensor,  as  discussed  in  Section  4. 
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4.  COMPARISON  OF  SATELLITE  SENSORS 


The  relationship  between  upwelling  radiance  and  the  atmospheric 

(3) 

aerosol  content  was  originally  demonstrated  using  Landsat  data.'  '  Since 
then,  this  study  has  found  similar  relationships  for  sensors  on  NOAA-5, 
GOES-1  and  SMS-2. ^  However,  the  sensitivities  of  the  systems  are  quite 
different.  NOAA-5,  which  is  no  longer  operating,  was  not  suitable  for 
aerosol  monitoring  due  to  a  non-random  noise  problem,  and  the  GOES/SMS 
sensors  are  only  6-bit  systems  and  do  not  provide  useful  radiance  sensiti¬ 
vity  for  monitoring  small  aerosol  content  changes.  The  most  useful 
satellites  currently  in  orbit  are  NOAA-6  (10-bit  AVHRR),  and  NIMBUS  7 
which  has  the  8-bit  CZCS,  designed  specifically  for  measuring  the  small 
radiances  over  oceans.  Unfortunately  digital  data  for  the  CZCS  are  not 
routinely  recorded,  and  are  not  readily  available,  whereas  the  NOAA-6 
(to  be  replaced  by  NOAA-7  in  late  1980)  data  have  been  routinely  archived 
by  NOAA  since  early  1980.  Thus  the  NOAA-6  AVHRR  is  the  most  useful 
operational  satellite  for  accurately  monitoring  atmospheric  aerosols  on  a 
routine  basis. 
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5.  INFLUENCE  OF  AEROSOL  PROPERTIES  ON  THE  RADIANCE-AEROSOL 
CONTENT  RELATIONSHIP 


(4) 

Theoretical  calculations  with  the  Dave'  '  atmospheric  scattering 
code  have  been  used  to  determine  the  effect  of  aerosol  properties  on  the 
radiance-aerosol  content  relationship.  These  properties  include  the  aerosol 
size  distribution,  real  and  imaginary  refractive  indices,  and  the  vertical 
distribution  of  the  aerosols.  The  calculations  have  been  made,  assuming 
that  the  aerosol  size  distribution  is  represented  by  the  Junge  distribution: 

dn(r)  =  Cr"v  d  log  r 

where  n(r)  is  the  number  of  particles  with  radius  r,  C  is  a  constant  depend¬ 
ing  on  the  number  of  particles  per  unit  volume,  and  v  ranges  between  2.5  and 
4  for  typical  size  distributions. 

5.1  LANDSAT  DATA 

Calculations  of  the  effect  of  v  on  the  radiance-aerosol  content 

relationship  have  been  made  for  each  of  the  Landsat  MSS  wavelengths  for  the 
(51 

standard'  '  Landsat  viewing  conditions  of  a  nadir  look  angle  and  a  sun 
zenith  angle  of  63°.  The  results  (for  a  refractive  index  of  1.5),  given  in 
Figure  5,  show  that  the  influence  of  v  varies  with  wavelength  and  has  a  mini¬ 
mum  effect  near  0.8  pm.  Thus,  the  use  of  a  satellite  sensor  at  this  wave¬ 
length  should  give  the  best  estimate  of  the  aerosol  content.  It  was  found 
in  the  Landsat  study ^  that  indeed  the  MSS  6  (0.75  pm)  showed  the  best 
radiance-aerosol  content  relationship  with  the  least  scatter  of  data  points. 
The  spectral  variation  of  the  influence  of  v  on  the  radiance  also  suggests 
that  spectral  measurements  of  the  upwelling  radiance  could  provide  information 
on  the  size  distribution;  wavelengths  greater  than  about  0.6  uni  would  be  pre¬ 
ferred  for  such  an  analysis  in  order  to  minimize  subsurface  effects.  The 
CZCS  is  preferred  to  the  MSS  for  this  analysis  since  it  can  measure  the  ocean 
radiances  more  accurately,  and  thus  more  readily  detect  the  small  radiance 
changes  due  to  changes  in  the  size  distribution. 
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Radiance  (mw/cnr/ym/sr) 


Wavelength  (uni) 


Figure  5.  Theoretical  Radiance  Versus  Wavelength  as  Function  of 
Aerosol  Content  and  Size  Distribution  (Nadir- viewing 
with  Sun  Zenith  Angle  of  63°). 


It  was  shown  previously^  for  the  standard  Landsat  viewing  con¬ 
ditions  that  increasing  the  real  part  of  the  refractive  index  increases  the 
upwelling  radiance  for  a  given  aerosol  content,  and  that  increasing  the 
imaginary  component  reduces  the  radiance.  In  addition  it  was  found  that  the 
vertical  distribution  of  the  aerosols  did  not  significantly  influence  the 
radiance. 

5.2  OTHER  SATELLITE  DATA 

The  effects  of  the  aerosol  properties  described  above  were  deter¬ 
mined  for  Landsat  MSS  observations  in  the  nadir.  However,  for  the  other 
satellites  discussed  in  Section  3  and  4,  the  sensors  scan  across  the  orbital 
path  and  experience  a  variety  of  sun  angles.  Thus  the  viewing  and  sun  angles 
differ  significantly  from  those  for  the  Landsat  observations,  resulting  in 
different  scattering  angles;  hence  the  aerosol  properties  may  cause  different 
effects.  To  investigate  this  possibility  a  computer  code  developed  by  us  for 
NOAA^  to  analyze  NOAA-6  data  has  been  utilized. 

The  NOAA-6  code  enables  the  aerosol  content  (0-10N)  to  be  determined 
from  a  radiance  measurement  by  the  AVHRR  0.65  pm  channel,  given  the  viewing 
angle  and  the  sun  angle.  The  code  is  based  on  radiance  calculations  with  the 
Dave  code  using  an  aerosol  model  which  best  fit  the  Landsat  MSS  6  (0.65  pm) 
data  at  San  Diego,  viz.  refractive  index,  n  =  1.5,  v  =  3.5,  and  an  ocean  albedo 
of  0.015.  The  code  considers  a  range  of  sun  zenith  angles  (eQ  =  42  -  84°), 
suti  azimuth  angles  (I  =  140  -  180°),  and  viewing  angles  (0=0-  84°)  to  cover 
all  conditions  to  be  found  in  the  N0AA-6  orbit  when  scanning  away  from  the  sun; 
these  values  may  not  cover  all  situations  for  other  satellites.  For  this 
investigation  the  calculations  have  been  repeated  for  three  different  sets  of 
aerosol  properties: 

Set  1.  n  =  1.33,  v  =  3.5 

Set  2.  n  =  1.50-0.01i,  v  =  3.5 

Set  3.  n  =  1.50,  v  =  2.5 

and  compared  with  the  standard  set  of  radiance  values  in  the  NOAA-6  code,  to 
obtain  a  tabulation  of  the  percentage  radiance  difference,  as  a  function  of  sun 
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angles,  viewing  angle  and  aerosol  content,  for  each  set  of  aerosol  properties. 

The  values  of  n  represent  the  extreme  values  which  are  typically  found  in  the 

atmosphere.  The  values  of  v  were  chosen  to  determine  the  effect  of  a  change  of 

f5l 

1.0  in  v;  the  previous  Landsat  study'  '  showed  that  the  radiance  change  at  a 
given  aerosol  content  is  approximately  proportional  to  the  change  in  v,  in  the 
nadir  viewing  case.  Selected  results  of  the  comparisons  were  given  in  Table  1, 
showing  the  percentage  difference  in  radiance  from  the  Landsat  model  computed 
for  each  set  of  aerosol  properties  given  above. 

It  is  found  for  Set  1  that  the  change  to  the  lower  refractive  index 
produces  a  radiance  change  always  within  about  ±1-3%  of  the  percentage  change 
observed  at  the  Landsat  conditions  (6  =  63°,  6  =  0°)  for  a  given  aerosol 

content.  The  difference  increases  with  N  up  to  6N  and  then  is  approximately 
constant  up  to  ION.  The  percentage  changes  are  not  very  dependent  on  the 
angles  e  ,  o,  and  $>,  except  for  large  sun  and  viewing  angles  (>72°)  where  the 
changes  actually  decrease. 

For  Set  2  it  is  found  that  the  introduction  of  some  absorption  by 
the  aerosol  produces  changes  up  to  two  times  larger  than  found  in  Set  1,  and 
which  increase  as  N  increases.  The  changes  are  generally  about  1-3%  less  than 
the  percentage  change  observed  for  the  Landsat  conditions,  and  again  are  approxi¬ 
mately  independent  of  angle  up  to  larger  values  of  eQ  or  e  (>72°)  where  the 
changes  are  smaller. 

For  Set  3,  the  reduction  of  v  produces  a  more  dramatic  effect  on 
the  radiance  change  in  terms  of  its  dependence  on  angle.  The  magnitude  of  the 
change  is  generally  less  than  in  Sets  2  and  3,  but  can  be  positive  or  negative, 
depending  mainly  on  the  angles  0  and  $,  with  a  lesser  dependence  on  e  .  The 
difference  increases  with  N  up  to  about  4N,  and  then  decreases  up  to  ION.  The 
magnitude  of  the  error  is  about  the  same  as  for  the  Landsat  conditions,  but 
may  be  positive  or  negative. 

It  appears  from  the  above  results  that  the  effects  of  variations  in 
the  aerosol  properties,  are  not  significantly  different  for  a  scanning  sensor 
than  for  the  Landsat  MSS,  except  for  variations  in  the  size  parameter  v  which 
can  produce  radiance  changes  in  the  scanning  sensor  with  the  opposite  sign 
of  those  for  the  Landsat  MSS,  but  still  of  about  the  same  magnitude. 
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These  radiance  changes  have  been  translated  into  errors  in  deter¬ 
mining  the  aerosol  content  from  the  NOAA-6  code  when  the  aerosol  properties 
are  different  from  the  Landsat  model.  These  errors  are  shown  in  Table  2  for 
the  radiance  changes  given  in  Table  1.  It  is  noted  that  the  same  radiance 
change  can  produce  different  errors  in  aerosol  content  for  different  sun  and 
viewing  angles,  since  the  slope  of  the  radiance-aerosol  content  relationship 
is  angle  dependent. 

Similar  calculations  of  radiance  change  and  aerosol  content  error 
were  made  for  a  change  in  the  vertical  distribution.  A  layer  of  aerosols 
with  a  number  density  seventy  five  times  greater  than  normal  was  inserted  at 
an  altitude  of  5  km.  The  calculations  showed  that  even  with  this  extreme 
variation,  which  perhaps  could  be  found  in  a  Saharan  dust  outflow  over  the 
Atlantic,  the  scanning  sensor  would  experience  about  the  same  radiance  change 
(~+5%)  as  the  Landsat  MSS;  the  corresponding  error  in  aerosol  content  is 
approximately  +10%. 
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Table  2.  Percentage  Change  in  Radiance  for  Changes  in  Aerosol 
Properties  (Aerosol  Content  =  IN) 
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Table  3.  Percentage  Errors  in  Aerosol  Content 
(Produced  by  Radiance  Changes  in  Table  1) 
for  Changes  in  Aerosol  Properties  (Aerosol  Content  =  IN) 
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6.  CONCLUSIONS 


The  calculations  show  that  errors,  due  to  variations  in  the  aerosol 
properties,  in  inferring  the  aerosol  content  from  a  satellite  radiance  meas¬ 
urement  are  similar  for  nadir  viewing  sensors,  such  as  the  Landsat  MSS,  and 
for  scanning  sensors,  such  as  the  NOAA-6  AVHRR.  In  comparing  the  calculated 
radiance  changes,  due  to  changes  in  the  aerosol  properties,  with  the  scatter 
of  data  in  the  measured  radiance-aerosol  content  relationships  found  for  the 
different  satellites,^  it  appears  that  variations  in  the  aerosol  properties 
must  be  relatively  small  over  the  oceans. 
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